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As of today it is widely believed that the theory of the strong interactions in the Standard Model is 
given by Quantum Chromodynamics (QCD). Hadrons are formed as bound states of different quark 
flavors, with the gluons mediating the strong interactions. Due to the strong coupling of QCD, a 
perturbative approach is not feasible at low energies, where the light hadron spectrum is observed 
in Nature. Simulations of QCD performed numerically on a finite lattice discretizing space and 
time offer the possibility to study the strong interactions non-perturbatively from a first-principle 
approach. In this review, I will try to summarize recent achievements and efforts made in lattice 
QCD simulations to extract the masses and decay constants of the light hadrons, i.e., hadrons built 
from the three lightest quark flavors — the up, down, and strange quarks. For a review on heavy 
quark physics on the lattice see the review by Christopher Aubin [[[]] given at the Lattice 2009 
conference in Beijing. See also the reviews by Yasumichi Aoki [0], Chulwoo Jung [St], Vittorio 
Lubicz [Q], and Ruth Van de Water [||] , which naturally have some overlap with the topics covered 
here and for the status as of last year's conference, see the review given by Karl Jansen at that time 

This review is structured as follows: first I will discuss the light pseudo-scalar meson sector, 
paying special attention to chiral perturbation theory, which is used to guide the extrapolations. 
Besides discussing the extrapolations for the pseudo-scalar meson masses and decay constants, I 
will also cover results obtained for the light quark masses and for some of the low-energy constants 
of chiral perturbation theory. In Section || I will broaden the discussion to cover the complete 
light hadron spectrum. Again, now with emphasis on the baryon masses, the techniques for the 
extrapolation towards the physical point will be reviewed. The status of studies focussing on excited 
states in the light hadron spectrum and related topics are summarized in the remainder before I give 
some concluding remarks. 



1. Masses and decay constants of the light pseudo-scalar mesons 

To extract the masses and decay constants of the light pseudo-scalar mesons — the neutral 



and charged pions and kaons — (for the light vector-meson (p) see Sec. |2.2| ) from most current 
simulations with dynamical quarks, an extrapolation to the physical pion mass point still has to 
be performed. Nowadays a typical value for the lightest dynamical meson mass is 250 MeV, but 
first results exist for dynamical meson masses in the range of 160 to 190 MeV. The inclusion 
of dynamical fermions has become standard due to improved algorithms developed over the last 
few years and increased computing power available, leaving the quenched approximation, with its 
undetermined systematic error, obsolete for almost all investigations of lattice QCD. In most lattice 
QCD simulations either two or three dynamical quark flavors are included, although first studies 
with four dynamical fermions are being pursued as well. The two lightest quarks, the up and the 
down quark, are usually taken as mass-degenerate, while the third quark flavor, the strange quark, 
is assigned a heavier mass. For that reason it is common to speak of either Nf = 2 or Nf = 2 + 1 
dynamical quark flavor simulations. Here I will focus mainly on aspects related with the extraction 
of quantities at the physical values of the quark masses. In the parameter space of quark masses, 
this point is usually defined by the constraint that the masses of the pion and kaon, m % and mx, 
(the latter only if the strange quark is considered, i.e. in simulations with Nf ■ = 2 + 1 flavors) take 
their experimentally observed masses at that point. To define this point, besides m n and niK a third 
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quantity is needed to set the overall scale in the computation. Here some groups utilize the mass 
of the Q. baryon, because of its expected weak dependence on the up and down quark masses. 
Other choices include the pion decay constant and the scale parameters ro (Sommer scale) or r\ 
extracted from the static quark potential. These different choices also provide opportunities for 
cross-checks, which are needed since none of the above choices is believed to be unproblematic 
in current simulations. For example, the Sommer scale r is not known a priori from experiment, 
and different estimates from lattice simulations vary by as much as 10 per cent. In the case of the 
pion decay constant, the extrapolation to the physical point as discussed below may introduce some 
unwanted systematic uncertainty. Recently, the Q. baryon mass has been advertised to overcome 
these difficulties, since it is a experimentally well determined quantity with only mild dependence 
on the up and down quark masses, being made from three strange quarks. But it remains to be 
seen, how big finite volume effects are in the case of this heaviest baryon within the light baryon 
spectrum. 

Table [I] gives an overview of recent dynamical simulations, a detailed discussion will follow in 
Sec. |L2|. As one can see, several fermion formulations are used, which differ in their approach to 
the continuum limit and chiral properties. Wilson fermions and their improved versions are cheap 
to simulate but introduce additional chiral symmetry breaking due to lattice artifacts, twisted mass 
fermions can be seen as a special variant of improved Wilson fermions. The (improved) staggered 
fermion formulations are very cheap to simulate and leave a remnant of the chiral symmetry unbro- 
ken at the expense of introducing additional taste degrees of freedom, which have to be accounted 
for. A more improved chiral behavior at the expense of additional simulation cost is offered by 
the domain wall formulation for fermions or the overlap formulation. The latter even offers exact 
chiral symmeUy at finite lattice spacing but is also most computationally demanding. I will not go 
into further technical details of the ensemble generation or a cost comparison between these sim- 
ulations. See the review of current dynamical simulations given by C. Jung at this conference [Q] 
and references therein. All except one analysis is currently relying on chiral perturbation theory to 
perform the extrapolation from their simulated meson masses to the physical pion and kaon mass: 
the PACS-CS collaboration in their recent work followed a different approach, namely to reweight 
their ensemble generated with a lightest meson mass of approximately 160 MeV to the physical 
pion mass [0]. 



1.1 Chiral perturbation theory for the meson sector 

Chiral perturbation theory (#PT) is an effective theory to describe the spontaneous and explicit 
breaking of chiral symmetry as an expansion in the masses and momenta of the light mesons 
H |9], 10], see, e.g., [11, 12] for recent reviews. The mesons act as the fundamental fields in this 



theory. Depending on whether one assumes chiral symmetry in the massless limit of two (up, down) 
or three (up, down, strange) quark flavors (usually referred to as the chiral limit), one formulates 
either SU(2) x SU(2) or SU(3) x SU(3) £PT (in the following I will for short just write SU(2) 
or SU(3) #PT). In addition to spontaneous breaking and explicit breaking (due to the non-zero 
quark masses) of the chiral symmetry from S\J(N)l x S\J(N)r to S\J(N)v, also lattice artifacts 
can introduce additional symmetry breaking effects. It is possible to address the latter in the #PT 
analyses. Therefore a distinction between continuum and lattice #PT for a specific formulation of 
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N f 


a [fm] 


m PS [MeV] 


N"LO 




m PS L 


ETMC 


2 twisted m. 


0.08, 0.07, 


250-600 


2 


compl. SU(2) 


>3.0 






A AC \ A 

U.U5 — > U 










JLQCD 


2 overlap 


0.12 


290-750 


2 


SU(2) 


> 2.8 




2+1 overlap 


0.10 


320-800 


2 


SU(2),SU(3) 


> 2.8 


PACS-CS 


2+1 iWilson 


0.09 


160^110 


1 


compl. SU(2) 


> 2.3 








reweight — > 








MILC 


2+1 staggered 


0.09,0.06, 


180-380 


1 


rS^PT SU(3) 


>4.0 






0.045 -+ 




2 


compl. SU(3) 










180-540 


3,4 


analytic 




Aubin et al. 


2+1 stagg/DWF 


0.12,0.09 


240-500 


1 


SU(3) MA^PT 


>4.0 






^0 




>3 


analytic 




RBC-UKQCD 


2+1 DWF 


0.11,0.09 


290-420 


1,2 


compl. SU(2) 


>4 



Table 1: Overview of chiral fits for the light meson masses and decay constants performed by various 
groups. The table lists the number of dynamical fermion flavors Nf included in the simulation and the 
fermion discretization used (two in the case of a mixed approach) and the lattice spacing(s). An entry a — > 
indicates that the continuum limit has been taken at some point in the analysis, mps specifies the range of the 
dynamical lightest meson masses (lower partially quenched meson masses may have been included in the 
analysis as well). The specific ^PT extrapolation applied and the order thereof is indicated by the exponent 
n in N"LO. To indicate the importance of finite volume corrections (which have been included in all cases), 



the lower bound for mp§L is provided, too. See Sec. 1 .2 for references and details 



the fields on the lattice has to be made. Besides the quark mass parameters l m U( i = (ra u + m<i)/2 
and m s , the effective Lagrangian of #PT also depends on several low-energy constants (LECs, 
sometimes referred to as Gasser-Leutwyler coefficients) which include effects of the heavier quark 
flavors as well as high energy QCD interactions. In lowest order two LECs appear: B and / for 
SU(2) or Bo and fo for SU(3). The latter being the decay constant in the chiral limit (here I choose a 
normalization, such that the physical f % ps 130MeV). In higher orders, additional constants appear, 
usually denoted as L, and K,. One has to keep in mind, that the LECs of SU(2) and SU(3) differ, 
since the former include the effects of the strange quark as well. 

Lattice simulations, besides relying on #PT for the extrapolation to the physical point, are also 
able to provide valuable information on #PT. In contrast to experimental measurements where the 
quark masses are necessarily fixed to their values in Nature, they can be freely varied in lattice QCD 
simulations. Therefore, lattice QCD should be able to test the predictions of #PT as functions of the 
quark masses and — if all lattice systematics are well understood and under control — make predic- 

'i will only deal with the case of two mass-degenerate light quarks, since this is current practice in lattice simula- 
tions, although ^PT formulae for quarks with non-degenerate masses are available in the cited literature as well. 
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tions about the convergence region and extract low-energy constants relevant for phenomenological 
models. 

A natural choice would be to use SU(3) #PT for 2+1 dynamical flavor simulations and SU(2) 
for those with only 2 flavors. However, it turned out that the convergence at next-to-leading order 
(NLO) in SU(2) #PT is much better than that in SU(3) at NLO. Especially, it is questionable, 
whether the meson masses close to or above the physical kaon mass can be successfully described 
by NLO SU(3) #PT. This has first been observed by the RBC-UKQCD Collaboration [|3|] and was 
later confirmed by other groups, e.g. Jl4]. When applying SU(2) #PT to data from Nf = 2+1 



simulations, the LECs are obtained at the simulated heavy quark mass. Practically, this is not a 
drawback, since the simulated heavier quark mass can be tuned to lie close enough to the physical 
strange quark mass, usually within 10 to 15 per cent is achieved. To account for the remaining 
mismatch in the heavy quark mass tuning, either simulations at different heavy quark masses and a 
subsequent interpolation have to be performed, or reweighting in the heavy quark mass might offer 
a promising remedy Jl5|]. 

Table [l] also lists the different variants of #PT, which have been used in the analyses of the 
various groups. Besides the continuum #PT there are variants which include certain effects due to 
the fermion discretization on the lattice like rooted Staggered #PT (rS^PT) and Mixed Action #PT 
(MA^PT). The former takes into account the rooting procedure necessary in the case of staggered 
fermions, while the latter is applicable when different fermion actions are used in the sea and 
valence sector. I did not explicitly mention the use of partially quenched #PT (PQ^PT) in the 
table. Partially quenching refers to the situation, when in addition propagators with quark masses 
different from the dynamical masses are calculated. Those partially quenched quarks are referred to 
as the valence sector in contrast to the sea sector of the dynamical (i.e. unquenched) quark masses. 

In addition to these #PT based extrapolations, there are also several ones which could be 
best described as ^PT-inspired. Currently, the JLQCD Collaboration uses fit formulae based on a 
resummation argument, where in NLO and beyond the LO squared meson mass 2Bm q gets replaced 
by the measured meson mass squared (rap s ) and likewise the decay constant in the chiral limit / 
by the measured one fp$. Since they also replace the scale /r in the chiral logarithms introduced 
by the regularization 2 with (4jt/ps), effectively higher order contributions are resummed in an ad 
hoc manner. In the table I labeled this fit ansatz as "^"^PT. For a detailed description of their fit 
functions and comparison to standard #PT see [|l6|]. 

Until recently, #PT fits were only performed including the complete terms up to next-to- 
leading order (NLO) and leaving out higher order terms. In the typical range of pion masses up 
to 400 MeV, those omitted terms are usually estimated to have a 3-5% effect. Therefore, if one 
aims for results with smaller systematic errors, one will ultimately have to include higher orders as 
well. One approach is to just add the analytic terms (which are multiplied by higher order LECs) 
to the fit formulae, ignoring the non-analytic (logarithmic) contributions. This has to be viewed 
as a (practical) phenomenological ansatz, but eventually it will be unsatisfactory if the predictions 
of #PT are to be tested by lattice simulations or one wants to measure accurately and in a well 
determined way LECs up to a given order. For this one needs the complete formulae for the meson 



2 For example, a typical term beyond LO of the form ^j^=p l°g —jjr 1 ls replaced by E, log £, , where the squared ratio 
of measured meson masses and decay constants S, = {T^bi at me simulated quark mass m q is used. 
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masses and decay constants in continuum (PQ)^PT up to that order. Bijnens et al. JT7| , |T8| , 19, 2C] 
published those up to NNLO and also provide FORTRAN code for those fit functions upon request. 
The complete NNLO in #PT is now being used in the analyses by ETMC, MILC, JLQCD, and 
RBC-UKQCD (either for SU(2) and/or SU(3)). By including the complete NNLO the number 
of LECs which need to be fitted is increased by four in the case of SU(2) (two from the 0(p 4 ) 
Lagrangian jSSt and two from the &(p 6 ) Lagrangian «Sfg) or in the case of SU(3) by ten (four from 
Jz?4 and six from Jz?6)- In the partially quenched versions, there are 13 (5+8) or 15 (5+10) new 
LECs for PQ-SU(2) or SU(3), respectively. (For comparison: up to NLO there are only four (two 
from Jz? 2 and two from Jg? 4 ) in SU(2) or six (2+4) in PQ-SU(2), SU(3), and PQ-SU(3).) All groups 
report (independently), that at the moment their available data itself is not sufficiently sensitive to 
determine those additional LECs. In order to get meaningful fit results, currently all groups add 
priors or constraints for the new LECs from phenomenological estimates, which are available for 
the LECs originating from the J^4 Lagrangian. For the remaining NNLO LECs, such estimates 
are not available. The sources for the phenomenological NLO LECs are the pion scalar radius, %% 
scattering, pion charge radius or the axial form factor in % — > Ivy in the case of SU(2) #PT. See 
[21 ] for more details and references. In the case of SU(3) #PT the best available estimates originate 
from combined phenomenological fits using the measured masses and decay constants of the pions 
and kaons plus data from decays, cf. [21] and references therein. In general one should be 
cautious using such estimates as priors or constraints when the goal is to test the applicability of 
#PT by comparing it to lattice data. All those estimates had to be extracted under the assumptions 
that (NLO or NNLO) #PT is sufficient to describe the experimental data available from the pion 
and kaon measurements, i.e., that #PT is valid at the pion and, more questionable, the kaon mass. 

A remarkable observation in these NNLO #PT fits to lattice data is that for the meson masses 
the NNLO contribution is almost of the same magnitude as the NLO contribution, possibly indi- 
cating a poor convergence of the (asymptotic) series. This is illustrated in the left panel of Fig. [I], 
which shows a NNLO fit performed by the MILC collaboration |@ |||. ETMC [||] and RBC- 
UKQCD [ p5| ] reported similar observations, when fitting their data with the complete NNLO. The 
anomalously small NLO contribution (being of the same size as the contribution from NNLO) has 
to be investigated further. Especially, it has to be excluded that this is influenced by the priors used 
in the fits. Given the fact that also a fit only using NLO describes the data well (with a bigger NLO 
contribution), one would expect NNLO to be a magnitude smaller. Comparing the extrapolations 
to either NLO or NNLO, it should be mentioned, that the extrapolated values are more or less con- 
sistent. Including NNLO terms therefore only shifts part of the NLO contribution to NNLO, when 
the latter is included. For the decay constants (shown in the right panel of Fig. [If), a better conver- 
gence is observed. But here one should also keep in mind, that the NLO already describes a 20 or 
more per cent correction at NLO between the chiral limit and the physical point. The effect seen 
for the decay constant can be best described by saying that the addition of NNLO "straightens" the 
fit curve, making it look more linear in the region where data is available. Ultimately, simulations 
at lower meson masses will have to show, whether this linear trend continues or at which point the 
expected chiral curvature will appear. 

For completeness, it should be mentioned that as an alternative to the chiral fit formulae, 
also an analytic expansion (polynomial expansion) in the quark masses or meson masses around 
a non-zero point could be used. (This has been discussed under the name "flavour expansion" by 
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Figure 1: SU(3) complete NNLO plus analytic N 3 LO terms jPT fit for the squared meson masses (left 
panel) and meson decay constants (right panel) from the MILC-Collaboration [ |26[ f22| ), plots courtesy of 
C. Bernard. 



L. Lellouch in his Lattice 2008 plenary contribution [27].) Results for a complete (meson masses 
and decay constants) analysis presented this year did not use an analytic expansion for their main 
results, which is why I will not discuss this point further in this context. 

1.2 #PT analyses of lattice data 

In the following, I will briefly summarize the main features of the chiral fits presented at this 
conference, before comparing the results of those fits in the remainder of this section. It should be 
mentioned, that all results discussed below either include finite volume effects in their analysis or 
their data has been corrected for these effects beforehand. 

The MILC Collaboration, using Nt ■ = 2 + 1 dynamical staggered fermions presented results 



from fits to SU(3) \ p.6[ |22[ ] and SU(2) [ |23| ] rS^PT showing good agreement between those results. 
To have better control in their SU(3) extrapolation, also "artificially light" strange quarks are con- 
sidered, which includes aNf = 3 simulation with three mass-degenerate light quarks as well. In 
a first step, they only fit their data for meson masses in the range of 180 to 380 MeV using com- 
plete NNLO (with priors for the LECs appearing only in NNLO terms). Since the inclusion of 
taste breaking effects is not available at complete NNLO, the root mean square averaged masses 
are used in NNLO, which is justified when taste breaking effects are negligible at that order due 
to a fine enough lattice spacing. In a second step, where the LECs obtained in the first step are 
now fixed, the fit range is enlarged to meson masses up to 550 MeV (to include the physical kaon 
mass) and analytic-only N 3 LO and N 4 LO terms are added, see Fig. |l| for an example fit. Currently 
two different lattice spacings are used in the first step (0.09 fm, 0.06 fm) of the analysis and in the 
second step simulations at an even finer lattice spacing (0.045 fm) are considered as well before 
the continuum limit is taken. 

The RBC- UKQCD Collaborations fit their data obtained with Nf = 2+1 domain wall fermions 



to SU(2) #PT up to NLO in a range of 290 to 420 MeV for the light pseudo-scalar meson [ 15]. The 
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light-strange meson (i.e. the kaon if the light quark mass approaches the physical value of m uc j) is 
described in SU(2) #PT as well. Data obtained at two different lattice spacings (0.09 fm, 0.1 2fm) 
have been used and a continuum extrapolation (assuming 0(a 2 ) scaling) is included in the analysis. 
This is a continuation of their previous work [13], where only the coarser of the two lattice spacings 
was available. The simulations are done at a fixed value for the strange quark, which is reweighted 
during the analysis to fine-tune to its physical value. To estimate the systematic effect of neglected 
higher orders in #PT, a comparison with results from an analytic expansion is performed. First 
results are also available from fits using complete NNLO #PT [25], but currently those are not 
used for the final quoted values. 

The PACS-CS Collaboration presented a detailed comparison of NLO SU(2) and SU(3) #PT 
fits to their Nf = 2 + 1 improved Wilson fermions data at a single lattice spacing at last year's con- 



ference [ ]28p and published those as well [|14||. Both, SU(2) and SU(3) #PT have been considered 
and a similar observation to that from RBC-UKQCD [13] about the better convergence of SU(2) 
#PT has been obtained. As already mentioned above, for their current analysis they pursue a dif- 
ferent strategy, namely to reweight their lightest simulated point at a meson mass of approximately 
160 MeV to the physical pion mass [[7j|. 

The JLQCD Collaboration now also has data available from Nf ■ = 2 + 1 dynamical overlap 
fermion simulations [ |29| ] in addition to their previous simulation at Nf = 2 [jig], each at a single 
lattice spacing of 0.10 fm and 0.12 fm, respectively. They perform the extrapolation from their 
mass range of 320 to 800 MeV to the physical point by using NNLO SU(2) and SU(3) #PT 
(see explanation above), since in fli6| ] they argued that this gives the best description of their data. 

The ETM Collaboration presented results from their simulations with 2 flavors of twisted mass 
fermions [24, |3(J. (Preliminary results for 2+1+1 flavors have been presented as well [^TJ].) For 
their meson masses in the range of 250 to 600 MeV they use complete SU(2) #PT up to NNLO 
(including priors for NNLO LECs). The continuum extrapolation is taken from three different 
lattice spacings (0.05, 0.07, and 0.08 fm). To obtain their final quoted values and the error estimate, 
they perform several different fits (varying fit ranges, using either NLO or NNLO,. . . ) and finally 
average the results weighted by the quality of the fit. The strange quark is quenched in their main 
analysis and the kaon decay constant and mass (or the thereby extracted mass of the strange quark) 
are obtained in a partially quenched set-up [32, 33]. 

In addition, the overview table also contains the mixed action results presented already at last 
year's conference by Aubin, Laiho, Van de Waterffi^ measured with domain wall fermions on 
configurations generated with Nf ■ = 2 + 1 dynamical staggered quarks by MILC. A recent update 
of their work has been presented elsewhere [35]. Currently two lattice spacing (0.12 fm, 0.09 
fm) are used for the continuum extrapolation. Complete SU(3) MA^PT is used up to NLO, but 
higher-order analytic-only terms had to be added to obtain reasonable fit results. 



1.3 Results 

1.3.1 Decay constants 

In Table ||, I compiled the results for the pion and kaon decay constants as well as the meson 
decay constant in the chiral limit of SU(2) and SU(3) #PT, / and fo, respectively. Two collabo- 
rations use the pion decay constant to set the lattice scale, so no prediction for f n is available in 
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f* [MeV] 


Sk [MeV] 


/[MeV] 


fo [MeV] 


ETMC 


input 


158.1(0.8)(2.0)(1.1) 


121.57(70) 




JLQCD (2) 


119.6(3.0)(1.0)(: 


h6.4\ 
-0.0 J 


111.7(3.5)(1.0)( 


+6.0 \ 
-0.07 


JLQCD (2+1) 


input 


157.3(5.5) 


121(14) 


79(20) 


RBC-UKQCD 


122.2(3.4)(7.3) 


149.7(3.8)(2.0) 


113.0(3.8)(6.8) 


93.5(7.3) 


PACS-CS 


134.0(4.3) 


159.4(3.1) 


126.4(4.7) 


118.5(9.0) 


MILC 


128.0(.3)(2.9) 


153.8(0.3)(3.9) 


122.8(.3)(.5) 


111.0(2.0)(4.1) 


Aubin et al. 


131.1(1.3)(2.2) 


156.3(1.3)(2.0) 






PDG 


130.4(.04)(.2) 


155.5(.2)(.8)(.2) 







Table 2: Results for the pion ( f K ) and kaon (Jk) decay constants from the chiral fits and their experimentally 
measured values []3rj|], see also footnote 3. Also given are the decay constants in the chiral limit of SU(2) 
(/) and SU(3) ^PT (/q). Results in italics mark preliminary results. In the case of the JLQCD collaboration 
results from their earlier 2 Jig] ] and recent 2+1 [29, i^] flavor simulations are shown. The RBC-UKQCD 



results are from [|l5| ^5[] except for /o, which is is from their earlier work [OJ. Other results taken from: 
ETMC [|§ Hi, PACS-CS [jTJ, MILC [g, Aubin et al. [||] 



those cases. With exception of JLQCD Nf = 2, all the quoted values agree within errors with the 



experimentally measured value [p6[], although it has to be noted that in some cases the combined 
statistical and systematic error is as big as 7 per cent. The same picture emerges for the kaon decay 
constant, although combined errors here are 4 per cent at most. 3 

Figure |2| shows the ratio of decay constants fa /fx compared to the experimentally observed 
value [j36|]. In addition, preliminary results from the BMW Collaboration [^] and older results 
from NPLQCD @ and HPQCD @] have been included as well. While the BMW Collabora- 
tion performed a direct extrapolation of the ratio (instead of separately extrapolating f n and fx), 
NPLQCD combined an SU(3) extrapolation of fx with the experimentally measured value of f n , 
and HPQCD used chiral expansions with priors. An updated average for the ratio of decay con- 
stants was given by V. Lubicz at this conference [Q]: fx/ fn = L196(l)(10), which is especially 
interesting for the determination of the CKM matrix element ratio |V u d|/|Ks|> see the contribution 



of R.S. Van de Water [|5|] and also [|41|] for more details and other implications in CKM physics. 

1.3.2 Low energy constants 

The results for the SU(2) LECs lj and I4 are shown in Fig. ||. Those are of interest in phe- 
nomenological applications, e.g., the pion-pion scattering length, see [fl3|]. The LECs are as usual 
defined at the scale of the physical (charged) pion mass. Within the quoted uncertainties, no dis- 
tinction can be made between LECs from Nf = 2+1 simulations, which include the effects of the 
strange quark, and those which only use Nf = 2 and therefore do not account for the strange-quark 
effects. The lattice simulations confirm the phenomenological estimates ^4, 21] and in the case 



of 1?, are also able to provide a value with smaller uncertainty. An interesting remark should be 



made at this point: the phenomenological estimate of I4 =4.4 ±0.2 [44, 21] together with the re 



3 It is also remarkable to note that the value quoted by PDG in 2008 [ pq l shifted by roughly 2.5<7 compared to the 
previously quoted value in 2006 [0. 
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Figure 2: The ratio of the kaon and pion de- 
cay constants fa/ fx from lattice simulations (for 
details see text) compared to the experimentally 
measured value (PDG '08) [Bq]. Also shown is 
the PDG '06 value [[37]], see remark in footnote 
3. Other values taken from: Lubicz [Q], MILC 
(Hi, Aubin et al. [§5[], RBC-UKQCD [|TJ & 
JLQCD ||| ||], PACS-CS [Q, BMW ||§], 
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Figure 3: The LECs It, (left panel) and I4 (right panel) in SU(2) jPT (defined at a scale m K ) from different 
lattice simulations, compared to phenomenological estimates [p[p4], cf. also J2l|]. The RBC-UKQCD result 
is from their earlier work Q, other values are taken from: MILC [ ^ |f], PACS-CS @, JLQCD [||| ||], 
ETMC [00. 



quirement that SU(2) %?T at NLO returns the experimentally measured value for f n constrains the 
decay constant in the SU(2) chiral limit to / w (121.5 ± 1.0)MeV. 

1.3.3 Quark masses 

Since #PT describes the dependence of the meson quantities on the quark masses, it allows 
one to extract the (light) quark masses once a reliable fit has been achieved. Usually the (exper- 
imentally observed) neutral pion mass is used to define the point of the average up/down quark 
mass m u d = {m u + m&)/2 and the kaon mass to define the strange quark mass. 4 The quark masses 



This is a somewhat simplified statement given the complexity of what is nowadays standard in ^PT fits to lattice 
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Figure 4: Left panel: The quark masses m u d and m s renormalized in the MS-scheme at jj, = 2GeV. Right 
panel: The quark mass ratio m s /m„<i. For details see text, values taken from HPQC D [Kq l, RBC-UKQCD 
" BJ, Aubin et al. [g, MILC ||], PACS-CS @, JLQCD [|| QCDSF [|JeTMC [f^ff]. 



which enter a lattice simulation are bare parameters defined in the lattice regularization scheme 
and depend on the fermion action and lattice scale. Therefore, the masses have to be renormalized, 
commonly the MS scheme at a renormalization scale of pL = 2 GeV is chosen. This transformation 
can either be performed perturbatively (up to some given order) or non-perturbatively by measur- 
ing the renormalization factors for the conversion to a regularization independent (RI) scheme for 
specific operators (taking into account possible operator mixing) directly on the lattice. In the lat- 
ter case, the conversion from the RI to the MS scheme still has to be performed perturbatively. 
For more details on the renormalization of quark masses and operators in lattice simulations, see 
the plenary contribution by Y. Aoki at this conference [ph. In the following, I will quote quark 
mass results in the MS(/X = 2 GeV) scheme and indicate whether non-perturbative or perturbative 
renormalization techniques have been used. 

Figure |] shows a compilation of both m ua and m s quark masses and their ratio obtained from 



the analyses discussed above. Also included are previous results from the PACS-CS [14], JLQCD 



(Nf = 2) qifj], and QCDSF [J4j] collaborations, as well as the HPQCD result presented at this con- 
ference [46]. The latter work used a different approach, namely to extract the mass of the strange 
quark from the strange/charm quark mass ratio m s /m c . Excluding the PACS-CS and preliminary 
JLQCD Nf = 2+1 points, which currently do not provide an estimate for their systematic uncer- 
tainty, the data might show a slight trend to higher quark masses observed in Nf = 2 simulations 
(leaving out the effect of a dynamical strange quark) compared to Nf = 2 + 1 simulations. But for 
a definite statement, the (mainly) systematic uncertainties have to be reduced further. 

The mass splitting between the up and down quarks can be estimated, e.g., by incorporating 
electro-magnetic effects by Dashen's theorem and the violation of the latter, and estimating the ratio 
m u /m u d from the observed mass difference between the neutral and charged kaons, see e.g. [p7|]. 
Results have been presented by the MILC Collaboration [22] and Aubin et al. [35], see Tab. ||for 

data. E.g., in a complete analysis, one in general needs three input parameters to fix m u( j, m s , and the lattice scale I /a. 
Commonly, the pion and kaon masses are used plus a third quantity like ro, r\, niQ or f n . In the case the third quantity 
depends on m u( j and/or m s (e.g. f % or hiq) a global fit procedure has to be carried out. But still the quark masses are 
mainly influenced by the input meson masses, so that the simplified statement is justified. 
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m u [MeV] 


m d [MeV] 


m u /m d 


MILC 1.96(0)(6)(10)(12) 
Aubinetal. 1.7(0)(2)(2)(1) 


4.53(1)(8)(23)(12) 
4.4(0)(2)(4)(1) 


0.432(1 )(9)(0)(39) 
0.39(1)(3)(0)(4) 



Table 3: The up and down quark masses m u , m& (MS-scheme, fi = 2GeV) and their ratio from [ [22] , [35[ | 
(preliminary results). 
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Figure 5: The mass spectrum 
of the light hadrons from dy- 
namical lattice QCD simulations 
(from §@§||g| 
|54]], see text for details), com- 
pared to the masses and widths 
(solid and dashed lines, resp.) 
observed in Nature |E}6||. 



a summary of their preliminary results. The RBC-UKQCD Collaboration presented preliminary 
results for the mass splittings at this conference, where the electro-magnetic effects have been 
included in the lattice measurements of Nf = 2 and Nf = 2+1 QCD via quenched QED [48]. 



2. Hadron spectrum from lattice QCD 

After limiting the discussion to the light pseudo-scalar sector in the previous section, I will 
now include light vector and scalar mesons and light baryons as well, turning to the complete 
light hadron spectrum. The mass spectrum of the light hadrons has been extracted from lattice 
simulations by various groups up to now, providing a good cross-check between the different lat- 
tice fermion formulations, experimental inputs and extrapolation methods used. In Fig. |5| I show 
a compilation of some recent results for the light meson masses and octet and decuplet baryon 
masses from lattice simulations compared with the masses (and widths) observed in Nature []3q]. 
Obviously, the lattice simulation results are consistent among each other and reproduce the exper- 
imentally measured values. As already discussed in the previous section, to fix the quark masses 
(m u d and m s or just m U( j if only non-strange hadrons are considered) and the overall lattice scale, 
one needs three (or two in the case of non-strange hadrons only) input parameters. In all the anal- 
yses discussed here, the pion and kaon masses have been used for the quark masses, while choices 
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for the third quantity may be the mass of the nucleon, the Q., or the S, the pion decay constant or 
the quark potential (j-q, r\). 

In the following, I will highlight some details of the analyses leading to the results shown 
in Fig. |5[ The Budapest-Marseille-Wuppertal (BMW) Collaboration [49] calculated the spectrum 



from their 2+ 1 flavor simulation using improved Wilson fermions at three different values for the 
lattice spacing and used a polynomial extrapolation from their simulated lightest meson masses 
in the range of 190 to 650 MeV down to the physical point. The Hadron Spectrum Collabora- 
tion (HSC) presented first results from a study performed to tune the mass of the strange quark 



in their simulations Q5fJ, |55J]. Currently, they have only one lattice spacing available using 2+1 
flavors of anisotropic clover fermions. Also a polynomial expansion is used to extrapolate from 
their lightest meson masses at 370 - 1520 MeV down to the physical pion mass. The PACS-CS 
Collaboration has data at a single value for the lattice spacing with 2+ 1 flavors of improved Wilson 
fermions available very close to the physical point, namely light meson masses ranging from 160 



to 700 MeV. In their initial study Q14| ] they used a (short) polynomial extrapolation to the physical 
point, whereas at this year's conference preliminary results obtained by directly reweighting to the 
physical pion mass have been presented [{7|]. The Lattice Hadron Physics Collaboration (LHPC) is 
using a mixed action approach to study the hadron spectrum. Domain wall valence propagators are 
calculated on the 2+ 1 dynamical flavors of improved staggered fermions at a single lattice spacing 
with lightest meson masses of 300 MeV and above [M]. They studied various different ansatze for 



the extrapolation towards the physical point, cf. also [£(]] and the discussion in Sec. [Q| . The MILC 
Collaboration has results available from 2+1 flavor simulations with improved staggered fermions, 
currently including three different values for the lattice spacings and a lightest meson with a mass 
of 180 MeV. They use, depending on the quantity, either chiral or polynomial extrapolations [p4|. 
The European Twisted Mass Collaboration (ETMC) calculated the nucleon and A baryon masses 
(i.e. only non-strange quantities) in their set-up of 2 flavor twisted mass fermions [51, |7|] from 
three lattice spacings at light meson masses down to 270 MeV using a chiral extrapolation. 

In addition to the above results, which are included in the summary plot (Fig. |J), the 
RBC-UKQCD Collaborations and the QCDSF-UKQCD Collaboration also published (prelimi- 
nary) results for the nucleon mass. The RBC-UKQCD Collaborations showed preliminary results 
for the nucleon mass extrapolation from two values of the lattice spacing with domain wall fermions 



[58]. The current work of the QCDSF-UKQCD Collaboration focuses on the splitting of the octet 



and decuplet baryon masses in the case of SU(3) symmetry breaking with 2+1 flavors of improved 



Wilson ("SLiNC") fermions [59] and the study of the p and A resonances with 2 flavors of clover 



fermions [6C], see Sec. 2.2, where the nucleon mass has been used to set the lattice scale. 



2.1 Extrapolations for baryon masses 

As for the meson sector, the formulae for the baryon mass extrapolation to the physical point 
can be based on chiral symmetry arguments. One approach is heavy baryon chiral perturbation 



theory (HB^PT) [51 , 62], where the effective fields again are either the pions in SU(2) HB^PT or 
the pions, kaons, and the v\ in SU(3) HB^PT. In general, the quark mass dependence of a baryon 
mass M bary on in HB^PT reads 

M _ M (0) , w(l) , .,(3/2) , 

^"baryon — ^baryon ~ r JW baryon ' Jw baryon "!"•••> 
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where every term shows a scaling with the quark mass m q according to 

''"baryon " l q ' 

Therefore, the expansion parameter here is e ~ m n ( K ^ / A x rather than ~ m 1 , K as i n meson 

#PT (but in both cases, the NLO term is of order m 2 n , K ^/A%). Alternatively, analytic expansions 
around the physical point or the chiral limit are used for the extrapolation, cf. also the review on this 
subject given by A. Walker-Loud at last year's conference [|56|]. In the following, I will review the 
current status of baryon mass extrapolations done by the groups mentioned above and also review 
work being done to study the SU(3)-breaking effects in baryon masses. 

SU(2) HB^PT has been studied with 2+1 dynamical fermion flavors by PACS-CS (improved 
Wilson fermions) [14] and the LHP Collaboration (mixed action: dynamical improved Wilson 
with valence domain wall fermions) [ 54 ] and with 2 dynamical fermion flavors by ETMC (Wilson 
twisted mass) |7|] and QCDSF (clover- Wilson fermions) [60]. While the latter two claim to 
observe a good agreement between their data and the predictions of SU(2) HB^PT, the PACS-CS 
Collaboration reports the theory to have a small convergence radius and an extrapolation which 
misses the physical point. LHPC finds that the fits describe their data adequately, but the extracted 
axial and nucleon-A couplings from those fits are inconsistent with phenomenological expectations. 
It should be mentioned, that the LHPC data has not been corrected for finite size effects, which 
seem to have an important impact in the other analyses. 

Fitting their measured lattice data to the predictions of SU(3) HB^PT has been pursued by 
PACS-CS and LHPC. Both LO and NLO fits are possible, but again discrepancies are observed 
between the octet and octet-decuplet axial couplings (commonly referred to as @, and ^) ob- 
tained from the fits and phenomenological models or direct lattice calculations of these couplings. 

Analytic expansions for the baryon masses were used in the analyses of BMW [||], PACS- 
CS [Q], and HSC [ 50"t 55] and also in part in MILC's analysis. Those seem to work fine, if the 
available data is close enough to the physical point as can be confirmed by simulations performed 
directly at the physical point (or reweighted to this point, as has been presented in ^7§). But one has 
to keep mind, that now the simulated lattice volume has to be large enough to exclude finite size 
effects since currently all methods to correct for such effects in a volume that is too small rely on 
(NLO)*PT. 

In a recent publication, Jenkins et al. [ 53 ] studied the effects of SU(3)-breaking via baryon 
mass relations, which they compared to actual data from lattice simulations. Using the 1 /N c ex- 
pansion, one can establish such mass relations which are expected to be fulfilled at @{N C ) leading 
to effects of the order of 1300 MeV, and subsequently or ff(l/N c ) relations leading to effects 
of the order of 430 MeV or 140 MeV, respectively. By this method, it is possible to study the 
validity of the baryon mass relations and in turn the applicability of HB^PT in a systematic way. 
Those mass relations turn out to be valid in the expected range, although the same mismatch for the 
couplings in HB^PT compared with phenomenological estimates is observed. The QCDSF Col- 
laboration in their current Nf = 2 + 1 dynamical clover- Wilson simulations also studies the effects 
of SU(3) breaking [p9|]. They simulate at different sets of two light and one heavier quark flavor 
with the constraint of keeping 2m/ + nih constant, including the point mi = m^. 

To conclude, the status of extrapolations using HB^PT is somewhat unsatisfactory at the mo- 
ment, at least if SU(3) HB^PT is considered. Further investigations, possibly including NNLO 
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fp [MeV] 


fp [MeV] 


ETMC [j7|] 


239(18) 


159(8) 


RBC-UKQCD [01 




143(6) 


Hashimoto et al. [73] 


210(15) 




QCDSF [0] 




168(3) 



Table 4: The vector meson decay constants f p and f p (renormalized in MS, ji = 2GeV). 



and artificially light strange quarks could reveal an answer for the convergence radius and there- 
fore if the strange quark mass lies within this region. With the data available at this time, short 
polynomial extrapolations or reweighting techniques seem to be more successful to extract the 
baryon spectrum. Also it should be mentioned 5 that relativistic (also called covariant) baryon ^PT 
(B^PT) []64], is an alternative approach based on chiral effective Lagrangians [ |66j |67| ]. The 
QCDSF Collaboration studied, e.g., the nucleon mass behavior within B^PT including finite vol- 



ume effects [|68|, |69j. Also the ETM Collaboration compared their SU(2) HB^PT fits mentioned 



above with fits based on the predictions from relativistic B#PT, stating good agreement between 



the two [51]. 



2.2 p vector meson mass and decay constant 

The lightest vector meson, the p, is unstable in Nature, since the decay into two pions is 
allowed. In lattice simulations, where the pion masses are light enough {2m n < m p ), the p therefore 
has to be treated as a resonance. Its mass can be calculated from the finite volume dependence of 



the phase shift of the %% resonance. The QCDSF [ ]60| ] and the ETM Collaborations [|70p reported 
on their projects to calculate m p by the above method. 

Results for the p decay constant have been obtained by several groups, see the overview in 
Tab. |j. Listed there are the decay constants from the coupling to the vector current and the tensor 
current, f p and f p , resp., in the MS renormalization scheme at 2 GeV. Experimentally observed 
are f p ~ 208MeV from T~ decay [Q and f p0 ~ 216(5) from p° ->■ e + e~ . 

3. Excited states 

The extraction of properties of baryon states beyond the ground state is more demanding for 
several reasons. First of all, the signal has to be obtained from the sub-leading exponentials in the 
fit to the correlator, having a much weaker statistical signal. Since in most cases a straight-forward 
multi-exponential fit would fail for that reason, one has to come up with a more sophisticated 



approach better suited to the problem at hand. The variational approach [|75J, [75]] seems to be 
most successful to solve this problem, for other approaches and a comparison see, e.g., and 
references therein. In the variational approach instead of fitting a single correlator one uses a whole 
matrix constructed from several correlators. These correlators need to have a sufficient overlap 
with states one intends to extract. Possibilities to construct several such correlators are, e.g., using 



I am thankful to Thomas R. Hemmert for pointing this out to me. 



15 



Light Hadron Masses and Decay Constants 



Enno E. Scholz 



Figure 6: Energy spectrum of excited nucleon states from Bulava et al. [[78]], obtained with Nf = 2 flavors 
of anisotropic Wilson fermions at light meson masses of 420 MeV (left panel) and 580 MeV (right panel). 
(Plots courtesy of J.M. Bulava et al.) 



different operators or different smearing prescriptions for the fields. Here I will focus on the mass 
spectrum of excited light hadron states only. 

Recently, J.M. Bulava et al. published results on the excited state nucleon spectrum [f7^] ob- 
tained from Nf = 2 dynamical anisotropic Wilson fermion simulations at two different light meson 
masses of roughly 420 and 580 MeV. They construct operators in the irreducible representation of 
the octahedral group corresponding to different spins and parities and using suitable displacements 
of the quark fields in the baryon operators. Figure |6| (from [78]) shows their result for the 7=1/2 
baryon spectrum at the two different light meson masses (the latter is indicated on the plots by the 
dashed line as well). Identifying their lowest state in the positive-parity channel as the nucleon, 
they find a cluster of negative-parity states at around 1.5-1.7 times the nucleon mass in accordance 
with the pattern of physical states observed by experiment. The higher positive-parity states lie at 
energies of 1.8 times the nucleon or above, leaving open the question whether or not the lowest 
state in this cluster will come down eventually when the lightest meson masses are lowered, so that 
it agrees with the Roper resonance at 1.53 times the nucleon mass. 

In previous studies, the Roper resonance Af*(1440) often was difficult to isolate and/or turned 
out to be found at too high a mass compared to the experimentally observed mass. S. Mahbub and 
collaborators in their recent work utilizing quenched FLIC fermions found a dependence of the 
extracted excited state mass on the smearing levels used and therefore suggested that one should 
extract the excitations based on the variational approach using different smearing levels in the cor- 
relator matrix [7^, 80, 31]. It is their conjecture that previous work actually reported a superposition 
of states rather than the Roper [ 79 ] . Figure [7] (from [|8T|]) shows a comparison of the nucleon ground 
state and the Roper from various recent lattice determinations. 

The meson spectrum for low and high spin excitations has been studied by T. Burch and 
collaborators [p^]. Based on Nf = 2 dynamical clover- Wilson fermion configurations from CP- 
PACS with a lightest meson mass of 500 MeV, they analyzed the ground and excited states for 
low spin (0,1) using a variational approach with several different quark sources. For high spin 
(2,3), only the ground states have been extracted. See Figure [8] (from [82]) for a summary of their 
results from two different lattice spacings (0.2 fm and 0.15 fm) extrapolated to the physical point. 
Given the current uncertainties plus systematics from the rather high meson masses used in the 
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Figure 7: The nucleon ground state and Roper resonance determined from different lattice simulations. 
(Plot courtesy of D. Leinweber, taken from [pT|].) 




Figure 8: Energy spectrum of low (0,1) spin (left panel) and high (2,3) spin (right panel) meson states from 
Burch et al. [p2||, obtained with Nf = 2 flavors of clover- Wilson fermions at light meson masses of 500 MeV 
and higher. (Plots courtesy of T. Burch et al.) 

extrapolation and potentially large Wilson fermion chiral symmetry breaking, at this point it might 
be too early to draw a definite conclusion. To overcome the limitations of chiral symmetry breaking 
induced by lattice artifacts, the excited meson spectrum is also examined using Nf = 2 flavors of 



chirally improved (CI) fermions Q83| , p4| ], where good signals have been obtained in the meson 
sector. The baryon spectrum from the CI fermion study currently turns out to give masses which 
are too high. This might be caused by the lattice volume being too small to produce reliable results 
for baryons. 

There is an open question, whether or not the lightest scalar mesons, the a, fc, and ao(980) 
are tetra quark (qqqq) states. S. Prelovsek et al. are currently investigating the 1 = 0, 1/2, 3/2, and 
2 tetra quark channels by using dynamical CI fermion configurations and also quenched overlap 
fermions [ |85| , |86| ], finding indications for a strong tetra quark component in the case of the a and 
K states. 



4. Concluding remarks 

This review covered recent and current efforts of the lattice QCD community to extract the 
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masses and decay constants of the light hadrons from numerical simulations following a non- 
perturbative, first-principle approach to QCD. In the light pseudo-scalar sector over the recent years 
steady progress in simulating lighter quark masses approaching the physical point has been made. 
Also the understanding of the advantages and limitations of the chiral extrapolation has gained a 
lot of insight from this progress and I tried to highlight the current status and discussion of chiral 
fits for the pseudo-scalar meson sector. For the future, the inclusion of the complete NNLO terms 
should be further pursued, both to improve the precision of the fits and to gain more insight into 
the convergence of the chiral expansion. This year the first promising attempts in this direction 
have been presented but all require additional phenomenological input, which eventually should be 
avoided. First results available close to or reweighted to the physical masses offer now the possibil- 
ity to test the values predicted by the extrapolations. The light hadron spectrum as summarized in 
Fig. || beautifully demonstrates the success of lattice QCD showing that the many different fermion 
discretizations used, combine into a consistent "big picture", although here some issues about the 
extrapolation methods used in the baryon sector need to be sorted out in the future. The study 
of excited states on the lattice also looks very promising. The tools for the extraction are well 
understood and first results with dynamical fermion simulations have been obtained. 
Acknowledgments I would like to thank the organizers of the "Lattice 2009" conference in Beijing for their 
invitation to give this plenary talk. Many colleagues helped me with this task by providing their material and 
preliminary results, promptly and patiently answering my questions and participating in fruitful discussions. 
I would like to thank Constantia Alexandrou, Yasumichi Aoki, Claude Bernard, Tom Blum, Peter Boyle, 
Tommy Burch, Norman H. Christ, Christine Davies, Chris Dawson, Xu Feng, Hidenori Fukaya, Steven Got- 
tlieb, Thomas R. Hemmert, Gregorio Herdoiza, Roger Horsley, Chulwoo Jung, Takeshi Kaneko, Chris Kelly, 
Yoshinobu Kuramashi, John Laiho, Christian B. Lang, Derek Leinweber, Laurent Lellouch, Huey-Wen Lin, 
Matthew Lightman, Vittorio Lubicz, Chris Maynard, Robert D. Mawhinney, Craig McNeile, Jun-Ichi Noaki, 
Shigemi Ohta, Haralambos Panagopoulos, Sasa Prelovsek, Siebren Reker, Christopher Sachrajda, Gerrit 
Schierholz, Amarjit Soni, Cecilia Tarantino, Carsten Urbach, and my colleagues in the RBC-UKQCD Col- 
laboration, the Fermilab Lattice Group, and at the University of Regensburg. For suggestions and a careful 
reading of the manuscript I am thankful to Norman H. Christ and Christopher Sachrajda. 

I gratefully acknowledge the financial support provided by Fermi National Accelerator Laboratory, 
operated by Fermi Research Alliance, LLC, under Contract No. DE-AC02-07CH1 1359 with the United 
States Department of Energy. 

References 

[1] C. Aubin, Lattice studies of hadrons with heavy flavors, PoS LAT2009 (2009) 007 
[ |arXiv: 0909 . 2686 [hep-lat ]| ]. 

[2] Y. Aoki, Non-perturbative renormalization in lattice QCD, PoS LAT2009 (2009) 012 
[ |arXiv: 1005 . 2339 [hep-lat ]| ]. 

[3] C. Jung, Status of dynamical ensemble generation, PoS LAT2009 (2009) 002 [arXi v : 1 1 . 9 4 1 

[hep-lat ] ]. 

[4] V. Lubicz, K aon physics from lattice QCD, PoS LAT2009 (2009) 013 [a rXiv :1004.3473 
[hep-lat ] ]. 

[5] R. S. Van de Water, The CKM matrix and flavor physics from lattice QCD, PoS LAT2009 (2009) 014 
[|arXiv: 0911 . 3127 [hep-lat]|]. 



18 



Light Hadron Masses and Decay Constants 



Enno E. Scholz 



[6] K. Jansen, Lattice QCD: a critical status report, PoS Lattice2008 (2008) 010[arXiv:0810.5634 
[hep-lat ] ]. 

[7] Y. Kuramashi, Physical point simulations in 2+1 flavor lattice QCD, PoS LAT2009 (2009) 1 10. 
[8] S. Weinberg, Phenomenological Lagrangians, Physica A96 (1979) 327. 

[9] J. Gasser and H. Leutwyler, Chiral perturbation theory to one loop, Ann. Phys. 158 (1984) 142. 

[10] J. Gasser and H. Leutwyler, Chiral perturbation theory: Expansions in the mass of the strange quark, 
Nucl. Phys. B250 (1985) 465. 

[11] J. Bijnens, Chiral Perturbation Theory Beyond One Loop, Prog. Part. Nucl. Phys. 58 (2007) 521-586 
[ |arXiv:hep-ph/0 604 043| l. 

[12] J. Bijnens, Chiral perturbation theory in the meson sector, PoS CD09 (2009) 031 



[arXiv: 0909 . 4635 [hep-ph] 



[13] RBC-UKQCD Collaboration, C. Allton et al, Physical Results from 2+1 Flavor Domain Wall QCD 
andSU(2) C hiral Perturbation Theory, Phys. Rev. D78 (2008) 1 14509 [ arXiv : 8 4 .0473 
[hep-lat ] ]. 

[14] PACS-CS Collaboration, S. Aoki et al, 2+1 Flavor Lattice QCD toward the Physical Point, Phys. 
Rev. D79 (2009) 034503 [ |arXiv: 0807 . 1661 [hep-lat ft 

[15] C. Kelly, P. A. Boyle and C. T. Sachrajda, Continuum results for light had rons from 2+1 flavor DWF 
ensembles, PoS LAT2009 (2009) 087 [ |arXiv: 0911 . 1309 [hep-lat]] ]. 

[16] JLQCD and TWQCD Collaboration, J. Noaki et al, Convergence of the chiral expansion in 



two-flavor lattice QCD, Phys. Rev. Lett. 101 (2008) 202004 JarXiv: 0806 . 0894 [hep-lat]Q . 

[17] J. Bijnens, N. Danielsson and T. A. Lahde, The pseudoscalar meson mass to two loops in three-flavor 
partially quenched chiPT, Phys. Rev. D70 (2004) 111503 JarXiv : hep-lat /04 60 1% 

[18] J. Bijnens and T. A. Lahde, Decay constants of pseudoscalar mesons to two loops in three-flavor 
partially quenched chiPT, Phys. Rev. D71 (2005) 094502 JarXiv : hep-lat /050 10 14[ ]. 

[19] J. Bijnens and T. A. Lahde, Masses and decay constants of pseudoscalar mesons to two loops in 
two-flavor partially quenched chiral perturbation theory, Phys. Rev. D72 (2005) 074502 
[ |arXiv: hep-lat /050 600"4] ]. 

[20] J. Bijnens, N. Danielsson and T. A. Lahde, Three-flavor partially quenched chiral perturbation theory 
at NNLO for meson masses and decay constants, Phys. Rev. D73 (2006) 074509 
JarXiv: hep-lat 70 602 003^ . 

[21] J. Bijnens, Quark Mass dependence at Two Loops for Meson Properties, PoS LAT2007 (2007) 004 
JarXiv: 0708 . 1377 [hep-lat]| ], 

[22] The MILC Collaboration, A. Bazavov et al., Results from the MILC collaboration 's SU(3) chiral 
perturbation theory analysis, PoS LAT2009 (2009) 079 JarXiv: 0910 . 3618 [hep-la~tT| |. 

[23] The MILC Collaboration, A. Bazavov et al., SU(2) chiral fits to light pseudoscalar masses and decay 
constants, PoS LAT2009 (2009) 077 [ [srXiv: 0911 . 0472 [hep-lat it 

[24] ETM Collaboration, P. Dimopoulos et al., Scaling and ChPT Description of Pions from Nf = 2 



twisted mass QCD, PoS LAT2009 (2009) 1 17 [arXiv: 0912 . 5198 [hep-lat; 



[25] RBC-UKQCD Collaboration, R. Mawhinney, NLO and NNLO chiral fits for 2+1 flavor DWF 
ensembles, PoS LAT2009 (2009) 081 JarXiv: 0910 . 3194 [hep-lat]|]. 



19 



Light Hadron Masses and Decay Constants 



Enno E. Scholz 



[26] 

[27] 

[28] 

[29] 

[30] 
[31] 

[32] 

[33] 
[34] 

[35] 
[36] 

[37] 

[38] 

[39] 

[40] 

[41] 
[42] 
[43] 

[44] 

[45] 

[46] 



The MILC Collaboration, A. Bazavov et al., MILC results for light pseudoscalars, PoS CD09 (2009) 



007 [arXiv:0910. 2966 [hep-ph] ]. 



L. Lellouch, Kaon physics: a lattice perspective, PoS LATTICE2008 (2009) 015 



flarXiv: 0902 . 4545 [hep-lat]| ]. 

Y. Kuramashi, PACS-CS results for 2+1 flavor lattice QCD simulation on and off the physical point, 



PoS LATTICE2008 (2008) 018 [arXiv: 0811 . 2630 [hep-lat]] 



JLQCD and TWQCD Collaboration, J. Noaki et ai, Chiral properties of light mesons with 



N f = 2+1 overlap fermions, PoS LAT2009 (2009) 096 flarXiv: 0910 . 5532 [hep-lat]| ]. 
ETM Collaboration in preparation. 

R. Baron et ai, First results ofETMC simulations with Nf=2+1+1 maximally twisted mass fermions, 



PoS LAT2009 (2009) 104 [arXiv: 911 . 52 4 4 [hep-lat] ] 



B. Blossier et al, Pseudoscalar decay constants ofkaon and D-mesons from Nf=2 twisted mass 



Lattice QCD, JHEP 07 (2009) 043 [arXiv: 0904 0954 [hep-lat] 



C. Tarantino private communication. 

C. Aubin, J. Laiho and R. S. Van de Water, Light pseudoscalar meson masses and decay constants 
from mixed a ction lattice QCD, PoS LATTICE2008 (2008) 105 [ arXiv :0810.4328 
[hep-lat ] ]. 

J. Laiho, "light quark results from a mixed lattice action ", talk given at Chiral Dynamics, Bern 2009. 

Particle Data Group Collaboration, C. Amsler et al., Review of particle physics, Phys. Lett. B667 
(2008) 1. 

Particle Data Group Collaboration, W. M. Yao et al, Review of particle physics, J. Phys. G33 
(2006) 1-1232. 

Budapest-Marseille- Wuppertal Collaboration, A. Ramos, FK/Fpi in full QCD, PoS LAT2009 



(2009) 259 [arXiv: 1002 . 1665 [hep-lat]] 



S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage, fx/fw in Full QCD with Domain Wall 



Valence Quarks, Phys. Rev. D75 (2007) 094501 [arXiv : hep-lat /0 60 602 3 [ 



HPQCD Collaboration, E. Follana, C. T. H. Davies, G. P. Lepage and J. Shigemitsu, High Precision 
determination of the n, K, D and D s decay constants from lattice QCD, Phys. Rev. Lett. 100 (2008) 



062002 [arXiv: 0706 . 1726 [hep-lat 



C. T. Sachrajda, Kaons on the lattice, PoS CD09 (2009) 032 [ |arXiv: 0911 . 1560 [hep-lat; 
J. Noaki private communication. 

H. Leutwyler, Model independent determination of the sigma pole, AIP Conf. Proc. 1030 (2008) 



46-55 [arXiv: 0804 . 3182 [hep-ph] 



G. Colangelo, J. Gasser and H. Leutwyler, pi pi scattering, Nucl. Phys. B603 (2001) 125-179 



[ |arXiv:hep-ph/01030 8 8| ]. 

M. Gtickeler et al, Simulating at realistic quark masses: Light quark masses, PoS LAT2006 (2006) 



160 [arXiv: hep-lat /0 610071] 



HPQCD Collaboration, C. McNeile et al., Towards precise relativistic b quarks on the lattice, PoS 



LAT2009 (2009) 116 [arXiv: 910 . 2 921 [hep-lat] ]. 



20 



Light Hadron Masses and Decay Constants 



Enno E. Scholz 



[47] MILC Collaboration, C. Aubin et al., Light pseudoscalar decay constants, quark masses, and low 
energy constants from three-flavor lattice QCD, Phys. Rev. D70 (2004) 1 14501 
[ |arXiv:hep-lat/04 070"28| l. 

[48] R. Zhou and S. Uno, Isospin breaking in 2+1 flavor QCD+QED, PoS LAT2009 (2009) 181 
[ |arXiv: 0911 . 1541 [hep-lat]| ]. 

[49] S. Diirr et al, Ab-initio Determination of Light Hadron Masses, Science 322 (2008) 1224-1227 
[ |arXiv: 0906 . 3599 [hep-lat]| ]. 

[50] Hadron Spectrum Collaboration, H.-W. Lin et al, First results from 2+1 dynamical quarkflavors on 
an anisotropic lattice: light-hadron spectroscopy and setting the strange-quark mass, Phys. Rev. D79 
(2009) 034502 flarXiv: 0810 . 3588 [hep-lat][ |. 

[51] European Twisted Mass Collaboration, C. Alexandrou et al, Light baryon masses with dynamical 
twisted mass fermions, Phys. Rev. D78 (2008) 014509 [ |arXiv: 0803 . 3190 [hep-latH - 

[52] A. Bazavov et al, Full nonperturbative QCD simulations with 2+1 flavors of improved staggered 



quarks, arXiv: 0903 . 3598 [hep-lat] 



[53] S. Gottlieb private communication. 

[54] A. Walker-Loud et al, Light hadron spectroscopy using domain wall valence quarks on an Asqtad 
sea, Phys. Rev. D79 (2009) 054502 JarXiv: 0806 . 4549 [hep-lat ]\ 

[55] S. Cohen et al, Excited-Nucleon Spectroscopy with 2+1 Fermion Flavors, PoS LAT2009 (2009) 1 12 
[ |arXiv: 0911 . 3373 [hep-lat ]| ]. 

[56] A. Walker-Loud, New lessons from the nucleon mass, lattice QCD and heavy baryon chiral 



perturbation theory, PoS LATTICE 2008 (2008) 005 [arXiv: 0810 . 0663 [hep-lat] 



[57] C. Alexandrou et al, Nucleon form factors with Nf=2 dynamical twisted mass fermions, PoS 



LAT2009 (2009) 145 [arXiv: 0910 . 3309 [hep-lat 



[58] RBC and UKQCD Collaboration, C. M. Maynard, Light Baryonsfrom 2+1 flavor DWF QCD, PoS 
LAT2009 (2009) 091 flarXiv: 1001 . 5203 [hep-lat ]| ]. 

[59] W. Bietenholz et al, Results from 2+1 flavours of SLiNC fermions, PoS LAT2009 (2009) 102 
[ |arXiv: 0910 . 2963 [hep-lat]| ]. 

[60] G. Schierholz, "masses of stable and unstable hadrons at physical quark masses", talk given at Lattice 
2009, Beijing 2009. 

[61] E. E. Jenkins and A. V. Manohar, Baryon chiral perturbation theory using a heavy fermion 
Lagrangian, Phys. Lett. B255 (1991) 558-562. 

[62] V. Bernard, N. Kaiser, J. Kambor and U. G. Meissner, Chiral structure of the nucleon, Nucl. Phys. 
B388 (1992) 315-345. 

[63] E. E. Jenkins, A. V. Manohar, J. W. Negele and A. Walker-Loud, A Lattice Test of\/N c Baryon Mass 
Relations, Phys. Rev. D81 (2010) 014502 [ |arXiv: 0907 . 0529 [hep-lat]] ]. 

[64] J. Gasser, M. E. Sainio and A. Svarc, Nucleons with Chiral Loops, Nucl. Phys. B307 (1988) 779. 

[65] T. Becher and H. Leutwyler, Baryon chiral perturbation theory in manifestly Lorentz invariant form, 
Eur. Phys. J. C9 (1999) 643-671 JarXiv : hep-ph/ 9 90 1 3 8 4| ] . 

[66] V. Bernard, T. R. Hemmert and U.-G. Meissner, Cutoff schemes in chiral perturbation theory and the 
quark mass expansion of the nucleon mass, Nucl. Phys. A732 (2004) 149-170 
[larXiv: hep-ph/ 03071151]. 



21 



Light Hadron Masses and Decay Constants 



Enno E. Scholz 



[67] M. Procura, T. R. Hemmert and W. Weise, Nucleon mass, sigma term and lattice QCD, Phys. Rev. 
D69 (2004) 034505 [ |arXiv: hep-lat /030 9020| ]. 

[68] QCDSF-UKQCD Collaboration, A. Ali Khan et al, The nucleon mass in N(f) = 2 lattice QCD: 
Finite size effects from chiral perturbation theory, Nucl. Phys. B689 (2004) 175-194 
[ [arXiv: hep-lat /0312o"3~o| ]. 

[69] M. Dorati, T. A. Gail and T. R. Hemmert, Chiral Perturbation Theory and the first moments of the 
Generalized Parton Distriputions in a Nucleon, Nucl. Phys. A798 (2008) 96-131 
[ |arXiv:nucl-th/0703073t . 

[70] X. Feng, K. Jansen and D. B. Renner, Scattering from finite size methods in lattice QCD, PoS 



LAT2009 (2009) 109 [arXiv: 0910 . 4871 [hep-lat 



[71] D. Becirevic, V. Lubicz, F. Mescia and C. Tarantino, Coupling of the light vector meson to the vector 
and to the tensor current, JHEP 05 (2003) 007 [ |arXiv : hep-lat /030 102 |. 

[72] ETM Collaboration, K. Jansen, C. McNeile, C. Michael and C. Urbach, Meson masses and decay 
constants fro m unquenched lattice QCD, Phys. Rev. D80 (2009) 054510 [ arXiv :0906.4720 
[hep-lat ] ]. 

[73] K. Hashimoto and T. Izubuchi, eta ' meson from two flavor dynamical domain wall fermions, Prog. 
Theor. Phys. 119 (2008) 599-641 [ [srXiv: 0803 . 0186 [hep-lat ]\ 

[74] M. Gockeler et al, Meson decay constants from N(f) = 2 clover fermions, PoS LAT2005 (2006) 063 
[ |arXiv: hep-lat 7050 91 9^ . 

[75] C. Michael, Adjoint Sources in Lattice Gauge Theory, Nucl. Phys. B259 (1985) 58. 

[76] M. Liischer and U. Wolff, How to calculate the elastic scattering matrix in two-dimensional quantum 
field theories by numerical simulation, Nucl. Phys. B339 (1990) 222-252. 

[77] H.-W. Lin an d S. D. Cohen, Lattice QCD Beyond Ground States, arXiv : 7 9 . 1 9 2 
[hep-lat ] . 

[78] J. M. Bulava et al, Excited State Nucleon Spectrum with Two Flavors of Dynamical Fermions, Phys. 
Rev. D79 (2009) 034505 [ |arXiv: 0901 . 0027 [hep-lat ]| ]. 

[79] M. S. Mahbub et al, Isolating Excited States of the Nucleon in Lattice QCD, Phys. Rev. D80 (2009) 



054507[arXiv:0905.3616 [hep-lat; 



[80] M. S. Mahbub et al, Isolating the Roper Resonance in Lattice QCD, Phys. Lett. B679 (2009) 
418-422 JarXiv: 0906 . 5433 [hep-latTj ). 

[81] M. S. Mahbub et al, Low-lying positive-parity excited states of the nucleon, PoS LAT2009 (2009) 
118 JarXiv: 0910 . 2789 [hep-latTfa - 

[82] T. Burch, C. Hagen, M. Hetzenegger and A. Schafer, Low and high spin mesons from Nf =2 
Clover-Wilson lattices, Phys. Rev. D79 (2009) 114503 [ [srXiv: 0903 . 2358 [hep-lat]] ]. 

[83] C. Gattringer et al, Hadron Spectroscopy with Dynamical Chirally Improved Fermions, Phys. Rev. 
D79 (2009) 054501 [ |arXiv: 0812 . 1681 [hep-lat ]| ], 

[84] G. Engel et a l, Excited hadrons in N f = 2 QCD, PoS LAT2009 (2009) 088 [ arXiv :0910.2802 
[hep-lat ] ]. 

[85] S. Prelovsek and D. Mohler, A lattice study of light scalar tetraquarks, Phys. Rev. D79 (2009) 014503 
[ |arXiv: 0810 . 1759 [hep-lat]| ], 

[86] S. Prelovsek et al, Spectroscopy of light tetraquark states, PoS LAT2009 (2009) 103 
[|arXiv: 0910 . 2749 [hep-lat]|]. 



22 



